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Abstract 

In this paper we construct entire solutions Ug to the Cahn-Hilliard equation 
—e^A(—e^Att -|- W'{u)) + W"{u){—£^l^u + W'{u)) = 0, under the volume con¬ 
straint /]g 3 (l — Ue)dx = 4\/27r^, whose nodal set approaches the Clifford Torus, 
that is the Torus with radii of ratio l/\/2 embedded in as e —)• 0. What is cru¬ 
cial is that the Clifford Torus is a Willmore hypersurface and it is non-degenerate, 
up to conformal transformations. The proof is based on the Lyapunov-Schmidt 
reduction and on careful geometric expansions of the laplacian. 

Keywords: Lyapunov-Schmidt reduction; Cahn-Hilliard equation; Willmore surface; 
Clifford Torus. 
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1 Introduction 

The Allen-Cahn equation 


— = u — u^, (1) 

arises in several physical contexts, such as the study of the stable conhgurations of 
two different huids conhned in a bounded container hi. If u{x) is the density of one of 
the two huids at a point x G hi and the energy per unit volume is given by a function 
W of u, it looks reasonable to obtain stable conhgurations by minimizing the energy 
functional 


E{u) = [ W{u)dx 

Jq 

among all distributions fulhlling the volume constraint 

/ udx = m. (2) 

Jn 

If, for instance, W{u) = (1 — and m G (—1,1), any piecewise constant function 
taking only the values ±1 and satisfying ([2]) is a minimizer, irrespectively of the shape 
of the interface. Therefore this model is unsatisfactory, since it is very far from the 
reasonable physical assumption that the interfaces are area minimizers, so one replaces 
the energy by 


E.(n) = ^ 

We can see that there is a competition between the potential energy, that forces u to be 
close to ±1, and the gradient term that penalizes the phase transition. By minimizing 
this functional, we are looking for the physical interfaces in which the phase transition 
can occur. 

The minimizers Us of are solutions to the Euler Lagrange equation, that is dni. 
In order to see if the interfaces are actually minimal surfaces, it is interesting to study 
the asymptotic behaviour of the level sets {u^ = c} as the parameter e —)■ 0. It is 
useful to exploit the variational structure of the problem. It was shown by Modica and 
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Mortola that the energy seen as a fnnctional on and extended to be +cxo 

when the integrand is not an fnnction, F—converges to the functional 


E{u) 


cPer^{{u = 1}) if M = ±1 a.e. in 
+00 otherwise in 


in the strong topology of (see [23]), where c > 0 is a suitable constant. 

Moreover, Modica showed that, if Ug are minimizers of E^ under the volume con¬ 
straint 



u^dx = m, 


for some m G (—1,1), then there exists a sequence 0 such that converges to 

some function u in (see proposition 3 of |22jh Furthermore, Theorem 1 of |22| 

asserts that u = ±1 a. e. in fl, and the set 


= {x e 12 : u{x) = 1} 


is actually a perimeter minimizer between all the subsets E <zVt satisfying the volume 
constraint 


^ |12| +m 
' ' 2 ■ 

Further results about the relation between the minimizers of Eg, and the minimizers 
of the perimeter can be found in [22] and in [7], where Choksi and Sternberg also 
described the relation between phase transition theory and the study of a certain kind 
of polymers. 

Conversely, it is an interesting problem to understand if any minimal hypersurface 
can be achieved as the limit of nodal sets of minimizers of the Ginzburg-Landau energy 
Ee. 

The hrst result in this direction is due to Kohn and Sternberg (see [IS])- They 
considered a smooth bounded domain 12 C and, as an interface, a disjoint union 
of segments k meeting the boundary dfl orthogonally. They dehned uq to be locally 
constant on 12\ Uj /j, taking the values ±1, and constructed a sequence of minimizers 
«£ converging to Uq in L^(12). 

In [2S], Pacard and Ritore proved a more general result, that holds true for a 
larger class of interfaces. They started from a minimal hypersurface S in a compact 
Riemannian manifold M and, under suitable assumptions, they showed that it can 
be achieved as the limit as e —)■ 0 of nodal sets (that is 0-level sets) of solutions 
of the rescaled Allen-Cahn equation ([T]). These solutions Ug were constructed with 
techniques such as hxed point theorems and the Lyapunov-Schmidt reduction, and are 
not necessarily minimizers. 

As regards the hypersurface S, they imposed some restrictions. They required it 
to be admissible, that is the nodal set of a smooth function / : M —)■ M. In the sequel, 
we will set 

M+(E) = {p G M : f{p) > 0} and M-(S) = {p G M : /(p) < 0}. 
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Moreover, S has to be non-degenerate. In order to explain the notion of non-degeneracy, 
let us give the variational characterization of minimal hyper surfaces. A hypersurface S 
in a compact Riemannian manifold M is said to be minimal if it is a minimizer for the 
area functional, whose critical points are characterized by the Euler equation if = 0, 
where H denotes the mean curvature of S. In the sequel, the mean curvature ii of a 
hypersurface S embedded in will always be 

H = k\ 

where the fc/s are the principal curvatures. 

The second variation of the area functional is given by 

A"(E)[0,V'] = j Lo(t){y)'ip{y)da{y), 

where the self-adjoint operator 

Lo0 = -As0 - |A|V 
is called the Jacobi operator of S and 

\A\^ = kl + --- + kl_, 

is the squared norm of its second fundamental form. By dehnition, a minimal hyper¬ 
surface S is said to be non-degenerate if its Jacobi operator 

Lo : C'2’“(S) ^ 

is an isomorphism. For an introduction to these topics, see also [9]. 

Moreover, the results in [26] hold even if the potential W{t) = (1 —is replaced 
by a more general double-well potential, that is a smooth function W such that 

{ W{t) >0 for any f, 

W{t) = 0 if and only if f = ±1, (3) 

> 0 . 

To sum up, they proved the following Theorem. 

Theorem 1 ([26]). Let W be as in Let E be an admissible non-degenerate minimal 
hypersurface in a compact Riemennian manifold M. Then there exists £o > 0 such that 
for any 0 < e < Eq there exists a solution to the rescaled Allen-Cahn equation 

— e'^Aus + W {us) = 0 

such that —)■ ±1 on compact subsets o/M^(S). 

Anyway, despite several results lead to think that, in some sense, the nodal sets 
of the solutions to the Allen-Cahn equation resemble minimal surfaces, there are also 
solutions for which the nodal set is far from being minimal. For instance, Agudelo, Del 
Pino and Wei constructed axially symmetric solutions u = u{\x |,a; 3 ) in such that 
the components of the nodal set, for |a: | large enough, look like a catenoid (see H). 
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The Lyapunov-Schmidt reduction was also applied to the non compact case, to 
construct entire solutions to the Allen-Cahn equation in R® that are monotone in 
one variable but not one-dimensional, since their nodal set resembles the Bombieri-De 
Giorgi-Giusti graph, that is a minimal graph over R® that is not affine (see 0,0). 

This solutions are related to a famous conjecture of De Giorgi, that asserts that, at 
least for N < 8, any entire bounded solution |m| < 1 to the Allen-Gahn equation 

— Am = u — 

satisfying djsfU > 0 in the whole R^ must be one-dimensional, that is it must depend 
just on one euclidean variable, in other words u{x) = u{< a,x >), for some unit 
vector a G S^~^. The result by Del Pino, Kowalczyk and Wei shows that de Giorgi’s 
conjecture is sharp about the upper bound on the dimension. Up to now it is known that 
the conjecture is true in dimension N = 2 (see |12],[TT]) and N = 3 (see 0,[II]). The 
conjecture is still open in dimension 4 < A < 8, although notable progress was made 
by Savin (see [28]), that proved that the conjecture is true in dimension 4 < A < 8 
under the reasonable assumption that, for any x G R^“^, 

lim u{x , Xm) = ±1, 

xjv —>-±00 

that yields that these solutions are minimizers of the energy 

E[n)=l + 

7 k 3 I 4 

We are interested here in analogues of these results for the Gahn-Hilliard equation 

- + W'{u)) + W'\u){-e^Au + W'{u)) = 0, (4) 

with W satisfying ([3]). Note that, as in the case of Allen-Gahn, we rescale the equation 
in order to treat T-convergence. If, for instance, we study the equation in a bounded 
domain G C R^, it is possible to see that it is the Euler equation of the functional 

w,(„) = if “ e n H^n) 

1-l-cxD otherwise in L^(r2). 

As in the case of the functionals related to the Allen-Gahn equation, some T—convergence 
results are known about We- More precisely, the asymptotic behaviour of We as e —)■ 0 
is related to the Willmore functional 

W(m) = c / 

JdEnn 

where U = {w = 1}, if m = ±1 a. e., defined when the interface dE is smooth enough. 

The nodal sets of the critical points m of W are called Willmore hypersurfaces. The 
Euler equation satisfied by this kind of hypersurfaces is 

- Aj:H = - 2HK, 

E 2 
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where H is the mean curvature and K is the Gauss curvature of S = dE. In the sequel, 
the Gauss curvature K of hypersurface S embedded in will always be 

K = ki... kN-i- 


An equivalent form of the Willmore equation is 

-Aj:H + ^H{H^-2\A\‘^)=0. (5) 

The Willmore functional arises naturally in general relativity, since it is related to the 
Hawking mass, that is 


mni^) 


Area(S) 

levr ^ 


1 

IGtt 


W(E)). 


Here S can be interpreted as the surface of a body whose mass has to be measured. 
Furthermore, this functional is also appears in biology, under the name of Helfrich 
energy^ and it is used to describe the behaviour of some lipid bilayer cell membranes. 
For further details and references, we suggest to see HHUHlls]. 


In |3] Bellettini and Paolini proved the F — limsup inequality for smooth Willmore 
hyper surfaces, while the F —lim inf inequality is much harder to prove. Up to now it has 
been proved in dimension N = 2, 3 by Roger and Schatzle in 123, and, independently, 
in dimension = 2, by Nagase and Tonegawa in j25]. The problem is still open in 
higher dimension, while it is known that the approximation does not hold, in general, 
for non smooth sets, even in dimension N = 2. 


In view of these F—convergence results that establish a link between the Gahn- 
Hilliard functional and the Willmore functional, it is interesting to see if also the above 
counter-part is true. In other words, we try to answer the following question; given a 
Willmore hypersurface S, is it possible to construct a sequence of solutions of the 
Gahn-Hilliard equation (jl]) whose nodal sets approach S as £ —)■ 0? In the paper, we 
show that this result holds true if, for instance, S is the standard Glifford Torus, that 
is the zero level set of the function 

f{x) = (^V2+ \Jxl + xpj +X 3 -I. ( 6 ) 

It has been recently proved in [TH] that the Glifford Torus is the unique minimizer of the 
Willmore energy (up to confromal transformations) among surfaces of genus greater or 
equal than 1 . 

It is interesting to see that it is possible to construct these solutions in such a way 
that they respect the symmetries of the Torus, that is the symmetry with respect to 
the XiX 2 -plane and with respect to any rotation that hxes the xa-axis. 

Theorem 2. Let W he an even double-well potential satisfying Let S he the 

Clifford Torus. Then there exists £0 such that for any 0 < e < £0 there exists a solution 
Ue to m satisfying the volume constraint 

I (1 — ufjdx = 4\/27r^, (7) 
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with Us —)■ ±1 and dkUs —)■ 0 uniformly on compact subsets for 1 < k < 4. More¬ 

over, Us{xi, X 2 , Xs) = Us{xi, X 2 , —X 3 ) and Us{x) = Us{Rx), for any x = {xi,X 2 , x^) G 
and for any rotation R G SO{3) such that R{0, 0,1) = (0, 0,1). 

In the statement of the Theorem, we denoted 

S’*" = {x G : f{x) > 0} and = {x G : /(x) < 0} 

This result is a fourth order analogue of Theorem [T] by Pacard and Ritore (see [26]). The 
proof is based on the Lyapunov-Schmidt reduction, that is we split equation (jl]) into 
a system of two equations. The auxiliary equation will be solved by using the spectral 
decomposition of the linearized Allen-Cahn operator and the bifurcation equation will 
be solved thanks to the nondegeneracy of the Clifford Torus, up to conformal maps. 
For a more detailed introduction to the techniques developed in the proof, see section 
2 . 

In order to explain what we mean by nondegeneracy, we go back to the varia¬ 
tional definition of Willmore hypersurface and we consider the second variation of the 
Willmore functional, that is 


>V"(S)[0,?/>] = j Locffjda, 

where Lq is the self-adjoint operator given by 

L^<t> = Llcf + V^H) + 2(AVs0, + (8) 

2H < A, VV > +0(2 < A, > +\V^H\^ + 2HtiA^). 

Here we have denoted by (•, •) the scalar product induced by the metric g on Ti, 
indeed, for instance (V0, Vif) = g^^Hicfj, and by <•, • > the trace of the product 
of two matrices, so for instance < H, >= and A^^ = g^^g^^A^i- It is 

possible to find the explicit computation of the first and the second variation of the 
Willmore functional W in [18], section 3. This is the analogue of the Jacobi operator 
in the case of minimal hypersurfaces. In view of a result by White [30], the Willmore 
functional is invariant under conformal transformations of the Euclidean space, that 
is homotheties, isometries and Mobius transformations, i.e. inversions with respect to 
spheres. On the other hand, by Corollary 2, page 34, of [21], we know that its second 
variation is positive definite on the orthogonal complement of the space of conformal 
transformations, hence the kernel of Lq exactly consists of these transformations. 

Remark 3. In view of the above discussion, Lq is injective if restricted to the space 
of functions with zero average and fulfilling the symmetries of the Torus, that is the 
symmetry with respect to the xiX 2 -plane and with respect to all rotations o/M^ that fix 
the X 3 axis. 

In fact, by considering just functions with zero average we exclude non trivial 
homothethies. This constraint is equivalent to prescribe the integral of 1 — that is 
to impose 


/ (1 — Us)dx = 4\/27r^ = 2|S 
Jm.3 


7 


where |S +|3 = 2\/27r^ is the volume of the interior of the Clifford Torus, that is its 
3-dimensional Lebesgue measure. In principle, a Lagrange multiplier Ag should appear 
in our equation: Anyway this will turn out to be 0 (see Section 7). By imposing rota¬ 
tional symmetry and symmetry with respect to the plane X 1 X 2 we exclude non trivial 
isometries and Mobius transformations. 

Acknowledgments The author is supported by the PRIN project Variational and 
perturbative aspects of nonlinear differential problems. The author is also particularly 
grateful to F. Mahmoudi, M. Del Pino, M. Kowalckyk and M. Saez for their kind 
hospitality and for their precious collaboration. 

2 Some useful facts in differential geometry 

For 0 < £ < 1, we dehne the rescaled Clifford Torus as := • C ^ 

other words, = {y G : f^{y) = 0}, where f^iy) '■= e~^f{ey) and / is dehned in 

dSD. 

For 0 < r < \/2 — 1 and 0 < e < 1, we dehne the tubular neighbourhood of width 
r/e of Se as 

Vr/e = {x G : dist{x, Eg) < r/e}. 

On this neighbourhood of Eg, we introduce a new system of coordinates, known as 
Fermi coordinates. First we dehne 

Zg : Eg X {-r/e, r/e) Vr/s 

by the relation 


Ze{y, z) = e^Py{zp{ey)) = y + zv{ey), (9) 

where ^{sy) is the outward-pointing unit normal to the original Torus E at ey, that 
coincides with the the outward-pointing unit normal to Eg at ?/, and exp^ is the expo¬ 
nential map of at y seen as a point of If r is small enough, that is 0 < r < \/2-l 
in the case of the Clihord Torus, Zg is a diheomorphism. In other words, Zg is a change 
of coordinates on lA/g, and the coordinates {y,z) = Zf^(x) are known as Fermi coor¬ 
dinates of the rescaled torus Eg, or stretched Fermi coordinates of the Torus. 

Remark 4. Any function u : lA/e —)■ ffi can be seen as a function of {y,z). More 
precisely, we can consider the composition u*{y,z) = u{Z^{y, z)). In the sequel, with a 
slight abuse of notation, we will write u = u{y, z). 

Let us hx a point Co ^ ^ and a parametrization onto a neighbourhood R C E of 
Co, that is a smooth function 


Y :U 

on an open set U GMf such that Y (Co) = Co, for some Co ^ U- Then, setting 74 = e~^U 
and 14 = 1 fhe function 


Ye-.U,^ Rg 
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given by Y^{y) := e~^Y{ey) is a parametrization of In the sequel, we will denote 
by y the points in and hy y = Y^{y) the points in V^. For any \z\ < r we consider 
the surface 


{y + zv{ey),y (10) 

On this surface, we consider the parametrization 

^^(y, z) := Y,{y) + zv{eYe{y)). (11) 

In particular, X := Xi is a parametrization of := Si 2 , the omothetic surface to S 
at distance It is known that the tangent vectors {diX^{y, z)}i=i ^2 constitute a basis 
of the tangent space Tyj^zu{ey)^£,z-, that will be referred to as the standard basis. We 
dehne the coefficients of the metric of Y^e,z at ?/ + zu{ey) as follows 

9 e,ij{y, z) :=< diXeij), djXeij) >= gij{ey, ez), (12) 

where <•, • > denotes the scalar product of and i, j = 1, 2. The Laplacian on 
is given by 

. - djW^^^ 9 e{y,z)gl^{y, z)di) = gl\y, z)dij + h\{y, z)di, (13) 

i/det 9 e[y,z) 

where 

K{y^^) ■= djgl^iy^z) + ^dj{\ogdetge{y,z))gl^{y,z) (14) 

and 0 := {g^^)ij are the elements of the inverse of the metric. These quantities are 
related to the ones of through the relations 

9 e{y:z) = g^\ey,ez), 

Khj^z) = eh\ey,ez), 

with g^^ := g^ and 6* := h\. We dehne the second fundamental form at y+zv{ey) G T^e^z 
to be the linear application of the tangent space Ty^zv{£y)'^£,z info itself that, in the 
standard basis {diX^ij, z)}i=i^ 2 i is represented by the matrix 

^e,ij{y, z) < div{ey), djX^{y, z) > . (15) 

We introduce the mean curvature Hsigii z) of at y + zp{ey) as follows 

He{y,z) := {Ae)\{y,z) = gl\y, z)Ae^ij{y, z). 

In other words 


Hs{y, z) = hAy, z) + hAy^ A, 

where ksAy-i z) are the principal curvatures of that is eigenvalues of the matrix 
g^^{y, z)A^{y, z). Therefore we can see that the metric gs,ij{y, z) and the matrix repre¬ 
senting the second fundamental A^^Ay^z) form depend on the parametrization, while 
this is not the case for H^{y,z). Setting, as above Aij := Ai^ij and H := Hi, we have 
Ae,ij{y,z) = eAij{ey,ez) and He{y,z) = eH{ey,ez). 
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Lemma 5. For a function u : —)■ M of class , the Laplacian in Fermi coordinates 

is given by 


Au{y, z) = As,, z) - eH{ey, ez)d^u{y, z) + d^^u{y, z). 


(16) 


For the notation, see Remark^ 

Proof. For any y & and \z\ < r/e, splits into the direct snni of the tangent 
space to Eg ,2 and the one dimensional snbspace generated by the nnit normal i'{sy), 
that is = Ty^zu{ey)^e,z + I®- The vectors {diX^iy, z),i'{ey)}i=i ^2 constitnte a basis 
of = Ty+ 2 i,(gy)M^. The metric in this basis is given by 


The inverse is 


Geiy.z) 


9e{.y,z) 0 

0 1 


Ge^{y.z) 


0 1 


(17) 


(18) 


Here 1 < /, J < 3 and 1 < j < 2. The laplacian on in the metric Gg is given by 

1 


Am = 


^det G^{y,z) 


^J{^/dsiG~(y^Gl\y,z)^I) = 


Gi-^{y, z)diju{y, z) + djGi'^{y, z)diu{y, z) + ^aj(logdet Ge{y,z))Gi'^{y,z)diu{y,z). 


iij 


ijt 


Thns 


Ge\y^ z)diju{y, z) = gl\y, z)diju{y, z) + d,,u{y, z) 
djG[\y, z)diu{y, z) = djgl^{y, z)diu{y, z) 

]^dj{\og det Geiy, z))Gl\y, z)diu{y, z) = 
]^dj{\ogdei ge{y, z))gl^{y, z)diu{y, z) + ^^^(logdet ^g( 2 /, z))d,,u{y, z). 
To conclnde, we point ont that 

^a;,(logdet^g(l/, 2 ;)) = -He{y,z) = -eH{ey,ez). 


□ 


Exploiting the Taylor expansion of H of the mean cnrvatnre of a given hypersnrface 
provided by Del Pino, Kowalczyk and Wei (see i), we have that 

H{ey,ez) = ^ ^ -tzk\ey) ^ Hj{ey) := Y,]=i Hi^y) (19) 


Here ki{ey) := fcg,j(|/,0) are the principal cnrvatnres of the Clifford Toms S at ey. 
Therefore the Taylor expansions of the first and the second derivatives of H are 


H,{ey,ez) = ^Hj+i{ey), 

H,,{ey,ez) = + ^){£zy-^Hj+ 2 {ey)- 
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In the sequel, we will set H{ey) := Hi{ey), \A{ey)\’^ := H 2 {ey) and tr74^(j/) := Hs{ey). 

Now we need the Taylor expansion in ez of For our purposes, it is enough 

to know the terms of order zero and one, while we also need the term of order two in 
the expansion of H. For this reason, we prefer not to expand the full Laplacian on 
In fact, an expansion up to order one would not be enough, because we cannot neglect 
the terms involving trA^, while an expansion up to order two would be a useless effort, 
in fact it would involve the terms of order two of ,,, that will always simplify in our 
forthcoming calculations. Before stating Next Lemma, we recall that 


A 




where 


/A ^ ^ A {^/detg,{y)gl^{y)di) = gl^{y)dij + bl{y)di, 

V^et ge{y) 


( 21 ) 


g7{y) ■= g7{y, o) = g'^iey, o) = g"^{ey) ( 22 ) 

Kiy) ■= Kiy^ 0 ) = ^b\ey, 0) = eb^ey). 

It is possible to hnd similar computations in 1211 , where Mahmoudi, Sanchez and 
Yao treat the more general case of a fc dimensional submanifold in an N dimensional 
manifold. 


Lemma 6. For a function m : lA/e M o/ class C"^, for any y G for any \z\ < rfe, 

As,,,m = As^m + ez{a7{ey)dij + eb\{ey)di) 
F{ezf{a7{ey)dij + eb\{ey)di) + a^^{ey, ez)dij + ebi{ey, ez)di, 


where 


af := b\ := 2djA^^ + 2F^^.A*^' - 

“2 := \dzz~g'^{ey,A),b\ := ]^d^A)\ey,Q), 

everything evaluated at ey, and the remainders satisfy \(f'^{ey, ez)\, \V{ey, ez)\ < 
for some constant c > 0 depending on S. 

Let 0,-0 : S —)■ M be functions. Let us set (fi := di(f). We recall that, by the 
properties of the covariant derivative, 

= duA^^ + V\iA^^ + ViiA^\ 

= 0P' - FF^fc, 

where everything is evaluated at ey. Moreover, by Codazzi’s equation, 'VjA^^ = 
g^^VkA^j, so in particular, 

=2(AV(^,V^) 

afffj,, + = 2kh^d/>,, - 2T’;,A^^iPk + (23) 

-(Vsi/, VeV') = 2 < a, > +(VsV’, Vsif), 

where we have set 


< A, V^^/J >:= A*^VjV’ = + Fj^fc. 
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Proof. By ffTT]) and ffT^ . we can see that 

9e,ijiy, z) = gij + ez{< diY, dju > + < djY, diU >) + (ezf < diU, dju > . 

In the proof, it is understood that the geometric quantities of S are evaluated at ey. 
In view of ffT^ with 2 ; = 0, we have 

dyy = -A%Y 

therefore 

96,ij{y, z) = 9 ij - ez{ 9 ikA’] + gjkA'l) + {ezf < diU, djU >= (24) 

9ij — 2ezAij + < diU, dju > . 

In order to expand the Laplacian, we need the expansion of the inverse of the metric. 
It is useful to write it as 


9 ^ — L + M, 

with Lij = 9 ij and M = —2ezAij + {ezy < diU, djU >. Equivalently, = L{I+L~^M), 
hence 

gf^ = (/ + L-^M)-^L-^ = {I- L-^M + 0 {{ezf))L-^ = - L-^MP-^ + 0 {{ezf), 

thus 

9 e{y^ z) = 9 '^ + 2ezA^^ + 0((ezf). 
where A^^ = g’^^g^^A^i- Moreover 

logdei ge{y,z) = logdet 5 (£(|/) + tr(L"^M) + 0{{ezY) = logdetfif^ - 2ezH + 0{{ezY), 
so, since \dj {log det g)A'-^ = V^-A'-^, 

As,,. = {g'^ + 2ezA^^)d,j + e{djg^^ + 2ezd,A^^)di 
+£(^cIj(logdet 5 f) — ezHj){g^^ + 2 ezA^^)di + 0 {{ezf‘) = 

As, + ezi^ 2 A^^di, + e{ 2 djA^^ + 2r^+ 0 {{ezf). 

□ 

As a consequence, we have the following expansion of the Laplacian 

A = - eH{ey, ez)d^ + As, + ez{af{ey)dij + eh\{ey)di) (25) 

+{ezf{ai{ey)dij + eh\{ey)di) + a^^{ey, ez)dij + eV (y, z)di. 

Although fl25p looks nice, we prefer to look for the expression of the Laplacian in a 
slightly different system of coordinates. We £x a function 0 : E —)• M whose L°°(S) 
is less than 1/4 and we introduce a new change of variables, that is we put 

t = z-(t){ey). (26) 
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The expression of the Laplacian will be more complicated than fl25|) . but more ap¬ 
propriate for our purposes. The reason is that we know the kernel of the operator 
— (As^ -|- dtt) + W (n*(f)), that is the one dimensional space generated by while 

we do not know exactly the kernel (if any) of —-|- dzz) + W (^*( 2 ; — (j){ey))). 
Given a function 


/ : Eg X M -)■ M 

of class C^, it is possible to dehne 

f : Eg X M — )■ M 

by setting f{y, t) := f{y,z — (f){ey)). A computation shows that 

ft(l/G) = fz{y,z- 0) 

Uy, t) = My, z-(j))- e(pif^{y, z-cfy 
Uj{y, t) = fij{y, z-(j))- e(j)ifzj{y, z - (f)) - ec^ifzjiy, z - (j)) 

+£‘^(j)ijfziy, z-(j))+ e(j)i(j)jf^^{y, z-(j)), 

where (j) and its derivatives are evaluated at ey, thus, in these coordinates, the expres¬ 
sion of the Laplacian of a function u dehned in lA/s of class is given by 

^ = dtt -\- Mdij -|- eUdi -|- D = 9ft -|- A^^ -I- D, (27) 

where the operator D is given by 

D := -eH{ey, e{t + 4>))dt - e‘^A^(j)dt - 2eg^^idtj + Vs0p9ff (28) 
+e{t + (f))[Mdij + £b\di - + b\Mdt - 2eaMidtj + e^aMi4^jdtt} 

+e'^{t + dij + eb\di - e^{ai(t)ij + b\(t)i)dt - 2eai(j)idtj + e‘^aMi(t>jdtt} 

+a^Mij + eUdi — + V'4>i)dt — 2ea'‘Midtj + e‘^a'‘Mi4>jdtt- 

Here we have set H{ey,s{t + (j))) := H{ey,ez), a'^^{ey,e{t + (j))) =lf^{ey,ez), V{ey,e{t + 
(j))) = b\ey,ez) and all the geometric quantities of E are evaluated at ey. 


3 Functional setting 


3.1 Functions on Eg 

As hrst we dehne, for 0 < a < 1, the space G^’"(E) as the set of functions 0 : E —)■ M 
that are k times differentiable and whose fc—th partial derivatives are Holder continuous 
with exponent a. We endow these spaces with the norms 




j=0 


-1- e" sup 
p^g 


|V^0(p) — V^0(g)| 

d(p,g)“ 


We note that these norms depend on e, since this is the right scaling in order to obtain 
our estimates. Moreover, in order to treat Lq, we dehne the spaces of functions that 
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respect the symmetries of the Torus, that is the symmetry with respect to the X 1 X 2 - 
plane and with respect to any rotation that keeps the xa-axis hxed. To be precise, we 
set T{xi,X 2 ,X 3 ) := (xi, X 2 ,- 0 : 3 ) and 


50 , 3 ( 3 ) :={i?e 50(3) :i?e3 = 63 }, 

where 63 = ( 0 , 0 , 1 ), and we dehne 

0^’“(S), := {0 e 0'=’“(S) : 0(C) = 0(TC) for any C G S, 
0(C) = 0(i?C) for any R G 50 , 3 ( 3 )}. 


By the symmetries of the Laplacian, the gradient and the geometric quantities of S, 
one can show that Lq preserves the symmetries of functions 0 G 0^’“(S)s, that is it 
maps 0 ^’"(S )3 into 0 °’“(S) 3 . 

We note that 50 , 3 ( 2 ) ~ 50(2), in the sense that any matrix R G 50 , 3 ( 3 ) has the 
form 


for some rotation of the xiX 2 -plane R G 50(2). 

Moreover, Remark [3] can be rephrased by saying that the operator 


C : 0^’"(S), X M ^ 0°’“(S)3 X M 


dehned by 



is injective. In fact, if £(0, A) = 0, multiplying by 0 and integrating over S we get 



and hence, since Lq is positive dehnite on 



(29) 


we conclude that 0 = 0, so A = 0. 

Being C also elliptic and self-adjoint with respect to the scalar product 


< (0, A), (0,/i) >:= / 0(C)0(C)c?fr(C) + A/i, 


s 


it is actually an isomorphism with bounded inverse. 
In the sequel ,we will often use the notation 


Rfc(l/4) := {0 G 0'=’“(S), : |0|c^,.(2) < 1/4}. 
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3.2 Exponentially decaying functions on 

For any 5 > 0 and for any x G we define 

ips{.x) =:C(|a;|) + (l-C(k|))e^l^l, 
where y : M —)■ M is a C°° cutoff function such that 


cw = 


1 for f < 1 
0 for f > 2. 


Moreover, we introduce the weighted spaces 

;= {u e ; ||h5||c^,.(M3) < CX)}, 

where us '■= u^fs and is the space of functions whose forth derivatives 

are Holder continuous with exponent a. We point out that functions belonging 
decay exponentially with rate 5, and the same is true for their derivatives. 

This spaces are endowed with the norm ||m|= |lu^lwhere 

k 

11^1 |c''.“(R3) := ||V%||oo + [V^uja- 

j=0 

In order to construct solutions to (jl]) that respect the symmetry of the Torus, we need 
to introduce the spaces of functions fulhlling these symmetries, that is 

C^’"(R^) := {u G C 5 ’“(R^) : u{Tx) = u{x) , u{Rx) = u{x) for any R G S' 03 ; 3 ( 3 )}. 

Remark 7. ITe note that, for instance, if u & C'|’“(R^), then Au G (^^’“(R^). In 
fact, by definition, any u G C|’“(R^) satisfies u{x) = Ut{x), where Ut{x) := u{Tx). 
Taking the Laplacian, we can see that Au{x) = Aut{x) = Au{Tx), and similarly, if 
R G SOx,i{3) and we set ur{x) = u{Rx), then Au{x) = Aur{x) = Au{Rx). 


3.3 Functions on x R 

First we will show existence and uniqueness of the heteroclinic solution to the ODE 
—+ W (n*) = 0. The result is known, but since the proof is quite short, we report 
it for completeness. 

Lemma 8. Let W be an even double well potential satisfying Then there exists a 
unique solution n* to the problem 

r-n;' + hF'(n*) = 0 

< ^*(0) = 0 (30) 

I n* —)■ ±1 as t ^ ±oo. 


and this solution is odd. 
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It is known that, if W{t) = |(1 — is the classical donble-well potential, then 
= tanh(t/\/2). 

Proof. Let n* be the nniqne solntion to the Canchy Problem 

-< + hP'(n*) = 0 
^^*( 0 ) = 0 

«:(o) = yW(o). 

Let (a, b) be its maximal interval of dehnition, with a < 0 < 
w{t) = —n*(—f) is still a solution to the same Cauchy Problem, 
so it is enough to study n* in the positive half line and a = —b. 
by and integrating we have 


b. Since the function 
n* is an odd function. 
Multiplying the ODE 


\(vy = W{v.) + c. (31) 

Evaluating at f = 0, it is possible to see that c = 0. As a consequence, > 0 in 
(0,6). In fact, if we assume by contradiction that there exists a hrst to such that 
n^(fo) = 0, then lP(n*(fo)) = 0, so in particular n*(fo) = 1, but, by the uniqueness 
Cauchy Theorem, this implies that n* = 1 in a neighbourhood of to, & contradiction. 
As a consequence, it is possible to dehne 

I := limn*(f). 

t—¥b 

By monotonicity, we know that Z > 0. Now we want to rule out the case I = oo. In 
fact, it this were true, we would have < 0 near 0 and > 0 near b, so there should 
exist > 0 such that v^{ti) = 0. Therefore, using the equation and (l3T]) . we can see 
that n*(fi) = 1 and = 0, which is not possible. 

Since Z < oo, we have b = oo. Now, always by dM]), we get that v'^ —)■ ^/2W{l) as 
t ^ oo. Since u is bounded, W{1) = 0, hence Z = 1. 

Uniqueness follows from the Cauchy Theorem. □ 


It is known that n* converges exponentially to ±1 as Z —)■ ±oo at a rate which is 
given by '^1U"(1) = ^1U"(—1), since W is even. More precisely, for any k eN, there 
exists a constant Ck such that 


for any Z > 0 

(32) 

and 



|af(n* + l)|<Cfce*V^^ 

for any Z < 0. 

(33) 

For instance, in the classical case W{t) = i(l — 

Z^)^, we have y/lU"(±l) 

= ^/2. 


For 0 < 5 < ^1U"(1), we dehne the function 

^s{t) = {l + e^Y{l + e-y. 
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For 0 < £ < 1 and 0 < a < 1, we define the space x M) as the set of fnnctions 

U : Eg X R —)■ M that are k times differentiable and whose fc—th partial derivatives are 
Holder continuous with exponent a. This space is endowed with the norm 


11^1 Icj’“(S£xK) ~ I |c''’“(SeXK)) 

where 


ll^llc'=.“(SeX]R) 


k 

^ ||V^t/||L-(SeXR) +sup 
j=o 


\x — 1 /|“ 


Given the heteroclinic solution n*, we can dehne the spaces 

X R) := |t/ e G 5 ’"(Se X R) : y U{y,t)v[(t)dt = 0 for any y e 

of functions that orthogonal, for any ?/ G Ee, to 

Moreover, as above, we will be interested in the spaces of functions that respect the 
symmetries of the Torus, thus we dehne 

(E, X R) := {U G G^“(E, X R) : 17t = G, Ur = U for any R G ^ 0 , 3 ( 8 )}, 

where we have set Uxiy, z) := U(Ty, z) and UR(y, z) := U{Ry, z). Furthermore, we set 
X R) := X R) n ^^’“(Eg x R). These spaces consist of functions that 

are both symmetric and orthogonal to v[. 

In the sequel, we will often mention the operator 

C,U := -(As, + du)U{yR) + W'\v,{t))U{yR), 

dehned for any U G ^^’"(Ee x R). 


4 Idea of the proof: Lyapunov-Schmidt reduction 

By a rescaling argument, it is enough to construct solutions to 

- A(-An + W'{u)) + W'\u){-Au + W’{u)) = 0, 

whose nodal set is close to Eg, since we can obtain the required solutions to (jl]) by 
setting u{x) := u{x/e). Thus we set 

F{u) = -A{-Au + W'{u)) + W"{u){-Au + W'{u)). (34) 

A computation shows that 

F'{u)v = -A{-Av + W"{u)v) + W"{u){-Av + W"{u)v) (35) 

+W'"{u){-Au + W'{u))v 

and 

F {u)[v,w] =—A{W {u)vw) + {W {u)W (u)+ W^‘^\u){—Au+ W {u)))vw + (36) 

W'” {u)[w{—Av + W” {u)v) + v{—Aw + W” {u)w)]. 
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In order to produce the required solutions we fix e > 0 small and a small function 
(j) G in the sense that |0|c4.“(s) < 1/4, and we define the approximate solution 

Ve^cj, in such a way that its nodal is exactly 

■= {y + (Pi£y)i^i£y ): y e SJ, 

and = ±1 outside a sufficiently small tubular neighbourhood of that is a 
neighbourhood of width r/2e + 6. More precisely, we set 


f 1 if fe{x) > 0 

H(a;) := < 0 if fe{x) = 0 

[-1 if /£( x )<0 


and, for any e: > 0 and for any integer m > 0, 


Xm{x) 


Cd^l - iix = Z^{y,t + (l){ey)) G K/e, 

0 if X G 


and we look for an approximate solution of the form 


VsA^) = X5(a:)Xe,</,(?/, f) + (1 - X5{x))Il{x), (37) 

where t is defined in fl26|) . and Ve,(j, is understood to coincide with El outside the support 
of X- Moreover Vs,<p will vanish exactly on and it will respect the symmetries of the 
Torus. We stress that these cutoff functions actually depend on 0, but we prefer not put 
the subscript 0 to simplify the notation. However, we will see that the error is 

small, but not zero, therefore we have to add a correction w = depending on e and 
0 in order to obtain a real solution, that is + ta) = 0. Rephrasing our problem 

in this way, the unknowns are 0 and w, for any £ > 0 small but fixed. Expanding F in 
Taylor series, our equation becomes 

F{Ve,cl>) + F'{Ve,cj,)w + Qe,<l,{w) = 0, (38) 


where 

Qe, 4 >{w) = dt F" {Ve, 4 > +sw)[w,w]ds, (39) 

Jo Jo 

However, we are not able to solve it directly, because of the lack of coercivity of F' 


4.1 The auxiliary equation: a gluing procedure 

We look for a solution of the form 

w{x) = X 2 {x)U{y,t)+ V{x), (40) 

where V is defined in the whole t/ is defined in the entire Eg x M. Since we want 
our solutions Mg to respect the symmetries of the Torus, we look for solutions U and V 
such that 


U{y,t) = U{Ry,t), for any R G and {y,t) G Eg x M 

V{x) = V{Rx), for any R G and x G M^. 


U{y,t) = U{Ty,t), 
V{x) = V{Tx), 
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Now we observe that the potential 


r,,^(x) := (1 - Xi{x))W"{v,,^) + Xi{x)W"{l) (41) 

is positive and bounded away from 0 in the whole that is, for any 0 < 7 < ^W”{1), 
0 < 7 ^ < <W (1) + To provided e is small enough, the estimate is uniform in 

<p. Moreover, using that X 2 X 1 = Xiy compute 

0 = X2!^F{ve,^) + + XiQsAU + V^) + XiM,,^(l/)|(42) 

+(-A + r,,^) V + (1 - x2)f{v,a + (1 - xi)QeAx2U + r) + N,,^(f/) + p,,^(i/), 
where 


M.AV) ■■= iW''iv,A - W'\l)){-AV + (43) 

+(_A + W"{v,A) - W"{l))V] 

A,AU) := -2 < Vx 2 ,V{-AU + W'\v,AU) > -Ax 2 {-AU+ W"{veAU) (44) 

+(-A + w"{v,A){-^ < VX2, VU > -AX 2 U) 
F,AV) ■■= -2 < Vxi,Vm"{v,A - W''{1))V) > -AxiiW'iveA - W"{l))v (45) 

+w"' + W AA)V. 

Hence we have reduced our problem to hnding a solution (V, U) to the system 

(-A + r,,^)V+(l-X2)i"(h.,0) (46) 

+ (1 - Xi)QeAX 2 U + H) + A,AU) + ^eAV) =0 in 

F{v,A + F\veA + XiQeAU + 4") + xMeAy) = 0 for |t| < r/2£ + 3. (47) 

The system of equations fl46|) and fITTj) is known as auxiliary equation. First we solve 
equation (146|l for any hxed U, thanks to coercivity, due to the fact that is bounded 
away from 0 uniformly in £ and 0. We will see that our solution also depends on the 
data U and e in a Lipschitz way. 

Proposition 9. For any £ > 0 small enough, for any U G CfA'^e x M) satisfying 
I |fo| lc^’“(SexR) — i- o.'iT'd for any 0 G i? 4 (l/ 4 ), equation admits a solution 14 ,</,,;/ G 
satisfying 

f \Ae,(t>,u\\cf°‘{R^) — 

\ ||14,0,Lri ~ fo,0,t^2llc^’“(i?3) A Cie “fo||17i — ^ 2 !|(74,a(s,,xR) (48) 

[||14,0i,c/ - V£,<p2,u\\cf‘^{R3) 4 Cie““fo|0i - 02|c-‘.“(s), 

for any Ui,U 2 satisfying ||foi|lc 0 “(SexR)> Il^alI c4.“(s,xR) ^ V any 0 i ,02 e ^ 4 ( 1 / 4 ), 
for some constants a,Ci> 0 independent ofU, e and 0. 

The proof of Proposition |9] is based on a hxed point argument (see section 6 ). 
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Now we consider equation fH7|) . In order to solve it, we need to extend it to the 
whole Eg X M. First we observe that 


F'{v,,^)U = ClU + R,,^{U), 


where 


R,,^(f/) := £,(D + W"{v,,^) - W''{vMU) + (D + W”- W”{v,))C,{U) 

+ (D + W'\ve, 4 >) - W'\v,)fU + + W\v,,^))U, 

D is dehned in ([28]). Therefore we reduced ourselves to consider 

ClU = -X4F(h.,0) - XiQsAU + V)- XaRsAU) - XiM.^y) (49) 


in the entire Eg x M. We would like to solve this equation with a hxed point argument, 
but, in order to do so, the right-hand side must be orthogonal to the Kernel of C^, that 
is the one dimensional space generated by n*(f), hence we can solve the problem 


cp 


-XiF{ve,^) - T(t/, 0) + p{y)v'^{t) 

poo 

/ U[y,t)v’^{t)dt = Q for any y e 


(50) 


where we have set, for the sake of simplicity. 


T(17, K, 0) := XiQeAU + V)- XaR^AU) - XiM,,^{V) 
piy) ■= — {xAF{ve,<i,) + T{U, Ve,^,u, 0 )) {v, t)v'^{t)dt 

J —oo 

and c* := 

Before stating the next proposition, let us observe that any function t/ : Eg xM —>■ M 
can be written as the sum of an even part and an odd part, the even part being 
Ue{y,t) := l{U{y,t)+U{y,-t)) and the odd part being Uo{y, t) := ^{U{y,t)+U{y,-t)). 

Proposition 10. For any £ > 0 small enough and for any 0 G i?4(l/4), we can find a 
solution t/g,</, G T^’"(Eg x M) to equation satisfying 

f I lc0“(SeXR) — ^2^^ 

\ ll(^£,</<)o|lc^’“(EeXR) — ^ 2 ^^ (51) 

[ - f4,</-2llc'0“(EeXR) ^ C'2^^101 - 02 |c4.“(E), 

for any 0 i ,02 G i?4(l/4), for some constant C 2 > 0 independent of e. 

The proof of Proposition [TOl will be given in section 6. 
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4.2 The bifurcation equation 

In conclusion, we will show that it is possible to hnd (j) that solves 

noo 

/ (\'4f (f'e,*) + T(f/, |^))(!/, t)li{t)dt = 0 


( 52 ) 


for any y E and such that the real solution Ue{x) := Ve,^{x/e) + We,(f){x/e) satisfies 
the volume constraint o. First we note that, by the change of variables x = x/e, 

47r^\/2 = / (1 - Ue{x))dx = e^ 1 - (v^^^ix) + w^^^{x))dx, 

Jr 3 Jrs 

the latter integral can be calculated exploiting the natural change of variables 


Xi = e ^ cos(£y 2 ) (( 2 ; + £ cos(£yi) + e ^\/2), 

X2 = sin(£y2) ((z + cos(£yi) + , 

X3 = e-'^{z + e:"^) sin(e:yi). 

on Vriei induced by the parametrization l^e(y) = £~^Y{ey), where 

Y ('di,'^ 2 ) := (cos'd2(cos'di + \/2),sin'd2(cos'di + \/2),sin'd2) 

and (di,'d2) =£(yi,y 2) e [0,27r)2. 


(53) 


(54) 


Proposition 11. For any £ > 0 small enough, c > 0 and 0 G (^^’"(S)^ satisfying 
|0|c4.“(s) < ce, 


1 - {vsA^) + ^ + 2e ^ / (j){C)da{C) 

Jt. 

i'6+tI2s 

/ t(l - n*(f))dt + 2Ge(0), 

'0 


+8V27iA-^ 


with Gs fulfilling 


|G,(0 )|<c, 

Aeifil) — ^£(02)1 < c|0i — 02 |c4.“(S)) 
for any 0, 0i,02 e C^’"(S)s satisfying | 0 |c 4 .c(s), | 0 i|c 4 .«(s), |02|c4,-(s) < ce. 

The proof of this Proposition will be given in Section 7. Therefore, in terms of 0, 
equation (j7]) is equivalent to equation 


0(C)dcr(C) = -4:V2Ae 


f6+T/2e 


f(l - v.(t))dt - 


(55) 


The system of equations (l52l) and fl55|) is known as bifurcation eguation, and it will 
be solved by a hxed point argument, that will be explained in this Proposition, whose 
proof will be carried out in Section 7. 


21 


Proposition 12. For any e > 0 small enough, the system of equations and 
/ f53]) admits a solution (j) G satisfying |0|c4.“(s) ^ C^e, for some constant 

Cs = Cs{W,t)>0. 

Remark 13. As we will see in the proof of PropositionlTR below, the Willmore equation 
will appear at order , while the linearized operator 

Locj) = Llf) + - R(Vs0, VsR) + 2(AVs0, VsR) + (56) 

2H < A, VV > +0(2 < A, > +1 VsRp + 2HtrA^). 

will appear at order thus it is crucial for the remainder to he smaller in order to 
apply a contraction mapping principle. This is actually the case thanks to the fact that 
the odd part of ^4,0 is of order 

5 The approximate solution 

5.1 Construction 

First one can try to take n*(t) as an approximate solution. We recall that t = 
z — 4>{sy), where 0 G i?4(l/4) is some small function that respects the symmetries 
of the S. We will see that these symmetries will be inherited by the approximate 
solution (see Remark [H] below). Since the Fermi coordinates are just defined in a 
neighbourhood of the Torus, our approximate solution is not defined everywhere. For 
our purposes, it is enough to consider it in the set 

B = {x = Z^{y, t + (j){ey)) G : |f| < tI 2e + 5}, (57) 

that is a tubular neighbourhood of 

^£,0 = {y + (j){ey)iy{ey) : ?/ G SJ 

of width r/Ae. Then it will be extended to the whole with the aid of a cutoff 
function. 

In the sequel, n* and its derivatives will always be evaluated at t, the geometric 
quantities, 0 and its derivatives will always be evaluated at ey. By fl25l) . 

— An* + kF'(n*) = —n* + kF'(n*) + eH{ey,e{t + 0))n* (58) 

+£^As0n* - £^|V0|^n* + e^{t + 0)(a*/0jj + b\(j)i)v[ - e^{t + (l))a'f (ficfjvl 

+e‘^{t + + 620 i)w - ^"^(^ + 0)^02^0j0iW 

+ V (j)i)v^ — cfiC^jV^. 

The term of order 0 in £ vanishes since n* satisfies the ODE —n” + kF'(n*) = 0. Thus, 
in order to compute P(n*), we need to apply the linear operator —A + W (n*) to the 
remaining terms. We will write down all terms of order less or equal than 4, the other 
ones being lower order terms, in some sense that will be clear soon. Let us set, for any 
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function v G C^(M), L*n := —v” + W''(Vi,)v. Differentiating the ODE satisfied by n*, 
we get = 0, thus using the Taylor expansion of H, the hrst term of (l58|) gives 

Tl^iy,t) = {-^ + W'{v^)){eH{ey,e{t + (p))v'j = e^{H^-2\A\‘^)v” (59) 

+e^!^{2H\A\^ - 4tryl3)(t + 0)^" + (^h\A\^ - 2tryl3)^;; _ Aj^Hv'^ 

+2(Vsi^, Vs0)i^r - i/|Vs0| V + 

+£^|(|A|'‘ -6H + 2HtiA^){{t + + (t + 0)nl) - As|Ap(f + cfyv'^ 

+2(vsiAp, VE0)(t+0)^: - iAnvE0r(t++As0iAr(f+ 0 )^: 

+ 0)nl + 2a^ Hi(l)j{t + 0)n" 

+//(a/0ij + 610i)(f + 0)n" - hfa/0i0jn"'| + f), 

with small and Lipschitzian in 0, in the sense that 

jl>S'e,<^^i,0lc0“(R3) < c 

- S^^^2^6a\c°-°‘{e) < c|0i - 02 |c4.“(s), 

for any 0, 0i, 02 G for some constant c = c(lT, r) > 0 independent of £ and 0. 

Similarly, the second term of fl58|) gives 

t) = {-A + w"{v^)) (£^As0nl) = e^HAj^(t)vl (61) 

+£^| - (As)+ |ApAs0(t + 0)n" + 2(VsAs0, Vs0)vr 

+(As0)V - |Vs0pAs0nr| +£®F2^(l/,t), 

with fulhlling (l60|) . 

The third term of fl55]) is already quadratic in 0, but, for the sake of completeness, 
we prefer to write it down. 

= ( - A + W"(v,))(-e^\^^ 4 'K) = - W"(v.)v:) ( 62 ) 

-E’'H\VEipfv“ +'5‘‘| - 1-41^1 VE.)>y( + + AsIVe^SIV 

- 2 (Ve:\Vj: 4 ,\\Vj:^)v: + - |Vj;.).rAj;.).K:"} 

The fouth term of fl5^ gives 

^) = ( - A + kT"(n*)) (£^(a/0p + &10i)(^ + 0)^^1) = (63) 

-2e^{a^(j)ij + b\(j)i)v” + e‘^H{al(t)ij + h\(t)i){v^ + {t + 0)n") + e^F^^iy, t). 

The hfth term of fl58|) gives 

Tl^hjA) = ( - A + lT"(n*))(-£^a/0i0j(f + 0)n") = 

£^a/0i0j((f + (t))v'f' -(t + 0)lT"(n*)n" + 2v”) 

-e'^Hal(t)i4>j{vl + {t + 0)n"') + e^F^^iy, f). 
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with fuiailing ea . 

Now we consider the terms involving and b^. We will see that all the contribu¬ 
tions of order coming from these terms will simplify, therefore we do not need to 
know the explicit expression of and 63 • 

^) = I ( - ^ + W'{v^)) {e^{di(t)ij + h\(t)i){t + 0)^^ - e^di^i(^j{t + (pfv”) (64) 

+e‘^{lf^(j)ij + V(j)i)vl — (pi(j)jvl'^ 

= -e'^{di(l)ij + h\(l)i){2v^ + 4(f 0)n") 

-e'^d^(j)i(l)j{2vl + 4(f 4- (j))v” + {t + -4 W”(n*)n") 

+ ( - A + W"(n*)) (^e^{-d^(j)ij + V(j)i)v'^ - f), 


with fulhlling (I6U|) . 

It turns out that, in the expansion of F(n*(f)), the only term of order is e^{H‘^ — 
2|Ap)n*. Since it is too large for our purposes, we add a correction to the approximate 
solution in order to cancel it. 

We set 

yit) = iv'Ss))-^ds 

This function is exponentially decaying, odd and solves 


LMt) 


-r] (t) + W {v^{t))r]{t) = -tv^{t) 


ri{t)v'^{t)dt = 0. 


Differentiating this relation once more, it is possible to see that Llri(t) = —n”(f). Our 
new approximate solution will be 


Ve{y, t) = v^{t) + e‘^{i){ey) -t- eL(t){ey))y{t), (65) 


with "0 : E —)• M and L linear in 0 to be determined later. In the sequel, rj and its 
derivatives are evaluated at t, the geometric quantities, 0 and its derivatives will be 
evaluated at ey. Taking the Taylor expansion of F^, 

F{veAy^ ^)) = F{v^) + F'{v^) + eL(j))y) 

+F"{v^) [e^{d + eL(l))y,e^{d + eLct))y\ + Ce,^[e'^{d{^y) + eL(t){ey))y\, 


where 

pX nt /*s 

Ce,</,[rc] = dt ds F (n* 4-rtc) [tc, tc, tcjdr. 

Jo Jo Jo 

Now we have to compute F'{vd){s^{d{^y) + ^F4>{ey))rj). As hrst we note that 


Tl 4 >{y, z) = IT"'(n*)(-An* 4- W'{vd))e^{d{^y) + ^L(t){ey))y = e^HdW'”{vj)yv^ 
+£^(0As0 + HL(j) +{t + J))d\A\‘^)W'”{v^)yv^ + e^Fl^{y, t), 
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with Fl^ fulfilling (|60|). 

After that, we have to compute ( — A + W” + eL(j))ri). We obtain 
( - A + + eL(j))r]) = e’^'ipL^r] + + LcjyL^r]) 

+£^| - Aj.'iPr] + + 0) + + 2(Vst/’, Vs0) 

+?/>As0)?7' -^/’|Vs0|V| + 

with satisfying flGOj) . 

Applying the operator once more, we obtain 

T^^ivA) = {- A + W''{v^)){e^{^p + eL(p)L^r]) = e^'il^Llr] + I^L(j)Llri + H^p{L^r])'^{66) 

+£^1 - AsV’^*h + + 0) + HL(j) + 2(VsV’, Vs0) + 'iPAy:4>) (L^v)' 

with satisfying fl60|) . 

Moreover, 

Tl^ivA) = { - A + W''{v,)){e^Hr]') = e^'iPHL,{7^') + (67) 

with satisfying fl6U]) . 

As regards the term of order of we note that 

( 2 /, t) = e%-A + W'\v,)) I - A^'iP'n + {\A\^'ip{t + 0) + FL0 + 2(VsV’, Ve 0)(68) 
+V^As0)V -V'lVs^lVj -'^sV’7^*h+ (FL0 + 2 (VsV’, Vs0) + 

+ |7lpV’^*((t + 0)V) - V'|Vs0pi^*(h")| 


with F^^^ satisfying fIBU]) . To conclude, also 

*) = (- A + «) 

is negligible, that is it satisfies (l60|) . since F^^ does. 

The only term of order in Fg(n *)+ eL(l)[ey))rj) is = —e^tljv^. 

Since we want it to erase the term of order of F^iVi,), we could set '0 := — 2|Ap. 

However, some quadratic terms appear at order e^. The only one that gives rise to some 
problems is —2F|Vs0Pw , thus we set 'ip := — 2|Ap + (i|Vs0P, for some constant 
d to be determined after projection. In particular, Vst/’ = 2FVsF — 2Vs|Ap + 
(iVs|Vs0P- L will be determined after projection. 
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Now we have to considered the contribution of F”(n*)(£^('0 + eL0)?7), since it gives 
rise to a term of order However, we will see that this contribution will cancel after 
projection 



with satisfying fIBU]) . 

We recall that is just dehned in B, while our global approximate solution is 
= X5(x46,^(y,t) + (1 - X 5 ( 2 :))IHI(x) (see ([37])). 

Remark 14. It follows from the construction that our approximate solution respects 
the symmetries of the Torus, that is '^£,44 = '^£,4^4) 

Re ^0,3(3). 

5.2 Projection 

As we noticed in section 4, 2, we need to consider the projection of the error F{ve,(j,)- 
In this subsection, we will explain how to do and we will see that this projection also 
enables us to choose L and d. 

Proposition 15. Let us set, for any (f e R4(r/4), 


Lcj):=-4:< A,V^(j)>+2HA^cj) + (j){2H\A\‘^-AtrA^, d=-ib4c,„ (70) 


where c* := J4oo4*4)Ydt. Then, for any y e Eg, the projec 

tion of F^Ve,4) satisfies 



(71) 


with uniformly hounded and Lipschitzian in (p ^ R4(r/4) and in e, that is there 
exists a constant c = c(W, r) > 0 such that 



(72) 


for any 0, e i?4(r/4) and for any £ > 0 small enough. 

Proof. Above we computed F^{ve,4 using fl27|) . now we just project it term by term. 
Integrating by parts we can show that 



(73) 


(74) 


'OD 


^ 0 , 


(75) 


— 00 
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so in particular 


W'”{v^{t))ri{t){v[{t)fdt = / [L^vit) - L^iv = -c*. 

./ _rv~\ ^ 


Moreover, setting 6* := {t)v'^{t)dt, we can see that 

/ OO 

- tW"{v^{t))v''{t) + 2v'''{t)}vl{t)dt = 

•OO 

poo poo 

- tLi,{vl{t))v[{t)dt-2bi, = - + 26* - 26* = 0 


( 76 ) 


because = 0. 

In the forthcoming calculations, right-hand side will always be evaluated at ey. By 
dZSI) and m, 

{Tl^iy^t) - e\H^ - 2\A\^)v:{t)]v:{t)dt = - c*As77 + 6*i7|Vs</-p} 

+£"c*| - 0As|7l|2 - (VslTlp, Vs0) - ^|7lpAs0 - 0(2 < A, V^H > +\V^H\^) 
-2{AV^H, Vs0) - Ih{2 < A, VV > +(Vsi7, Vs0)) | + e*" 


j Tl^iy,t)'v[{t)dt = £^1 - (As)V - ^iTlpAs^l + 

/ OO 

T!jy,t)v,(,t)dt = £='6.F|V!:.#.r + e^TI^, 

-OO 

/ OO 

^' 0 ( 1 /, t)v'^{t)dt = e^-c^H{2 < A, > +(Vsi7, Vs0)) + 

with satisfying ([72]). 

By ([teD, 

/ OO 

Tl^{yA)v\{t)dt = 

-OO 

/ OO 

- tW”Mt))vl{t) + 2v”{t)]v^{t)dt + 

-OO 

with 0 satisfying (172]) . Once again by (175]) . we can see that 

/ OO 

Tl^iyA)v:m = ^^JA,Et>, 

■OO 

with satisfying fl72]) . 
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Now let us consider the terms coming from the correction. 

- 2\A\-^) + dH\V^d>\'^\ 

+£-‘c.1|(//2 - 2|/1|2)As^ + HL<t, + (i/2 - 2|/1|2)|/1|2.^| + 

/°° {TUvA - - 'i\A?)L,m}v'Mdt = £2c.1|//(//2 - 2|/1|2) 

+dH\Vs<t>\^'^ + - 2\A\^)\A\^d, + 

-(Vd2l|2, V^4>} + - 2|A|2)As.#, j + 

with satisfying ([72]). To conclude, also 

/ OO 

{Tfjy,t) + T"‘^{y,t)]v',(t)dt (77) 

■OO 

fulhlls fl72|) . In conclusion, we choose L and d as in (ITOj) in order to cancel the quadratic 
term appearing at order and to obtain exactly Lq as a linear term at order Since S 
is a Willmore surface, that is it satishes the Euler equation —Aj^H + — 2\A\^) = 

0 , we have 

= c’lc.! - Ae// + 1//(//2 - 2|71|2))(£j,)| 

-£'^c*Lo0(e2/) + e^ne,4>iey) = 

-e‘^c,,Lo(j){ey) + e^'He,<p{ey), 

where 'He,<p '■= Y^t=i statement is true with '■= □ 

6 Solving the auxiliary equation 

This Section will be devoted to the proofs of Propositions jHl and [TUI In both cases, 
we will hrst study the linear problem associated to our equation and then we will apply 
a contraction mapping principle. 

6.1 Solvabilty far away from the linear problem 

We will prove the following Proposition. 

Proposition 16. Let 0 < 5 < 7 < ■^iy''(l). Then, for any £ > 0 small enough, for 
any 0 G i?4(r/4), and for any f G the equation 

(-A + P,,^) V = / (78) 

admits a unique solutionV = '^s,4 >{f) Cff {M.^) satisfying\\V\\^4,c < c| |/| 

for some constant c > 0 independent of e and 0 . 
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Remark 17. The symmetries of the solution follow for free from the symmetries of the 
laplacian and ofT^^^. In fact, if f G and V is a solution to (—A + = 

f, then also Ut{x) := u{Tx) is a solution, thus, by uniqueness, u = Ut- The same 
argument also shows that u = Ur, for any R G SO^^iS), hence u G 

We stress that the assumption 5 < 7 is crucial. When we solve the equation 
(—A + = / we lose some regularity, in the sense that the solution might decay 

slower than /. 

We split the proof into some lemmas and a proposition, with the aid of some re¬ 
marks. First we reduce ourselves to consider a second order PDF, then, by a bootstrap 

argument, we will solve our forth order equation. 

Proposition 18. Let 0 < 5 < 7 < ^JW''{!). Then, for any e > 0 small enough, for 
any 0 G i? 4 (r/ 4 ), and for any f G the equation 

- Am -h = / (79) 

admits a unique solution u = satisfying | |m| < c| |/| 

for some constant c > 0 independent of e and 0 . 

Before giving the proof, we state a technical Lemma. 

Lemma 19. For any 1 < p < 00 , there exists a constant C = C{p, 5) > Q such that, 
for any u G we have 


m||l 7 r 3 ) < C\\uips\\ oo* 


Proof. The case p = cxo is trivial, since (ps > 1, so we can assume that p < 00 . We split 
the norm of u into the sum of two terms, that is the integral over a ball of radius 
R> 1 and its complement 


\u\^dx = / \u\'^dx+ / \u\^dx. 


'Br 


'Bl 


The first term satisfies 


I \u\^dx < \Bji\ 

M • 

'Br 



p 

OO’ 


and the second one fulfills 

[ \u\Pdx=f {\u\ipsYip^_sdx < \\uips\\P^ f ip^_sdx <\\uips\\^^, 


iBi 


iBi 




for some suitable R = R{S,p) > 1, where we have set (fs = 1/(^5. 
Now we are ready to prove Proposition [TBl 


□ 


29 



Proof. Step (i): existence, uniqueness and local Holder regularity. 


Existence and uniqneness of the weak solntion follow from the Riesz representation 
theorem. Since / e then u e 

Step (a): estimate for the L°° norm. 

Now we will show that mps G and 

WuifsWoo < cWfcpyWoo. (80) 

From now on, we will assume that / is not identically 0, and hence Hm^Hoo > 0, 
otherwise there is nothing to prove. As hrst we will prove that ^ 0 as |x| ^ cxo. In 
order to do so, it is enough to show that uip.y G This will be done by using 

the function as a barrier. More precisely, we fix p > 0 and \z\ > p. Then we fix 

cr > 0 and R> \z\ so large that u{x) < cr for |x| > R. Therefore u fulhlls 

u < maxQBp u < < Xe~'^P + a for |x| = p 

^ u < a < Xe~'^^ + a for |a;| = R 

(—A + + a)) < <0 for p < |x| < R, 

provided A > Aq, with Aq independent of a. By the maximum principle we get that 

u(z) < + cr, for any \z\ > p and for any cr > 0. Letting cr —)■ 0, we get that 

G L°°. 

Since 0 as |a:| —>■ cx), the supremum is achieved at some point p G M^, that is 
llh^lloo = \u 5 {y)\. Now we consider two cases. If \y\ < 1, we apply the elliptic estimates 
to control the L°° norm of ups with f.y, otherwise we use the equation for us. 

Let us consider the case \y\ < 1. We observe that us = u in Bi{0) and we apply 
elliptic estimates to get that 

||w||l°°(Bi( 0)) < <^11^11^2.2(52(0)) < <^(lkl|L2(B2(0)) + II/IU2(B2(0))). 

Now we multiply the equation (—A + Tsp)u = f hy u and integrate by parts to obtain 
that 


||w||l2(B2(0)) a ||m||51(K3) < <^||/||l2(r3). 

Moreover, by Lemma [T^ applied with p = 2, we get 

||/||l 2(]R3( < C||/.y||(X). 

Exactly in the same way, we can estimate the term 11/| |l 2 (S 2 )- point out that all the 
constants are independent of £ and 0, because the potential is positive and bounded 
away from 0 uniformly in e and 0. 

Now let us turn to the case in which the maximum point of lu^l is achieved outside 
Bi. The equation satished by us is 

- Aus + = fs-"^ < Vm, Vps > -uAps = (81) 

fs - 2p_s < Vus, Vps > -us{2 < Vps, Vps > -psAps). 
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If Ui{y) > 0, then us has a global maximum point at y, and a computation shows that 

0 < 'j'^iisiy) < Te,<t,iy)u5{y) < -^U5{y) + I'^usiy) < fsiy) + 5‘^u{y), 

the estimate being uniform in e and in 0, hence 

My) < ^2 i - M(y)^ 

so we have (180]) . If us{y) < 0, then us has a global minimum at ?/, and the conclusion 
follows from a similar argument. 

Step (in): estimates for higher order derivatives. 

Now we will show that 


||'*^<5||c 2.“(R3) < c||/.^||c0.a(K3), (82) 

for some constant c > 0. In order to do so, we observe that, by elliptic estimates (see 
|13)h we have that, for any x G 

II“<5||c2.“(Bi(x)) < (||/(5||cO’“(B2(a:)) + I I'h.sl |L°°(B2(a;))) ^ c| |/.^| |c0.« (M^), 

the constants being independent of £ and 0. □ 

Now we can conclude the proof of Proposition [T6l 
Proof. Given / G G°’“(M^), we have to hnd V G G^’^(M^) fulhlling 


I (-A + TMM = f 
\ I 1^1 Ic^’“(R3) < C| i/I 1(^0, 

In order to do so, we use proposition fTSl twice to hnd w G and u G G^’/(M^), 

with 0 < 5 < 5' < 7 , such that 


(-A + r,,^)n = / 


and 


{ ll'*^llc^;“(R3) — *^1 I/I lc/“(R3) 

II^IIc/“(R3) — *^11"^! Ic°;“(R3)- 

Now it remains to estimate the higher order derivatives of u. For this purpose, we 
differentiate the equation satished by u and we get (—A + = Uj — for 

j = 1,..., 3, hence, applying the regularity estimates for (—A + Fg^^), 


I l^'l Ic/“(r 3) — *^(11'^jl Ic°;“(r3) + I I/I Ic°’“(r3)) 

< c(||nj<py ||co.<»(r3) + ||/|Ic°’“(R3)) a c(||'Uj</9y ||c1(r 3) + ||/|lc°'“(R3))) 
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hence 


||V (^</55)||oO < C(| |m| + I I/I |c0,a(-]g3j) < c| | / | | (^0,a ^]g3 j . 

<5 

Similarly, differentiating the equation once again, we see that 

+ ^e,(i>)Vij = Uij — 

for / j = 1,..., 3, so in particular 

I l^il lc/“(R3) — *^1 l'^*il lc°;“(R3) “ c(| |Mjj</9y I |c0.q(r3) + 11/| |( 70 .“(r 3 )) 

5 

< c(||-Uij</9y ||ci(R3) + ||/||^0.«(jg3j), 

therefore 

||V iy ips)\\oo + [fdc^^jo^Q, < c(| |m| |(^ 2 ,Q^jg 3 ^ + I I/I lc°’“(R3)) — *^1 I/I lc°’“(R3)! 

5 

all the constants being independent of £ and 0. □ 


6.2 The proof of Proposition [9]: solving equation (l46ll by a fixed 
point argument 

Equation fH6|) is equivalent to the fixed point problem 


y = Ti(E) := vl/,^^|(l - X2)Fyy + (1 - Xi)QeAX2U + E) + y^U) + P,,0(E)|, 

that we will solve by showing that Ti is a contraction on the ball 

Ai - {y e C'‘;,"(R=) : ||V'||c.,.,R3| < Cie-“''}, 

provided the constant Ci is large enough. In fact, by the exponential decay of U far 
from Eg, we get that 

l|N.,4f^)llc*-(R.) < 

for some constants a, c > 0 depending on W, r, 6 but not of £ and /. By fl32|) and fl33|) . 
the same is true for (1 — X2)F{vsA- Moreover, by fl58|) . fl3^ and fl33|) . 

||Pe,0(l^)|lc:>/“(R3) Fl ce ^ I |I/| |^4.q^jj 3^ < ce ^ , 

with c > 0 depending on IP, r, 5 but not of e and </. Moreover, using that 


“ Xi)l^llch“fR3) ^ c||P|U4,a 


and 


||(l-Xi)X2t^|lch 


< ce 


—aje 
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where (1 — Xi)X 2 U is understood to be 0 outside the support of X 2 , and the dehnition 
of (see ([39])), we get 

11(1 -< ce-^i‘. 

Up to now, we have just proved that Ti maps Ai in itself. In order to show that it 
is actually a contraction, we need to estimate its Lipschitz constant. The only terms 
depending on V are Pg,,/,, that fulhlls 

I ~ -P£,(/)(U2) I < C£||Vl — h^l 

for some constant c > 0 independent of £ and 0 , and {l — xi)Qe,<t>i.X 2 U + U), that fulhlls 
11(1 - xMQ.Ax^U + V)- < ce-‘l‘\\V, - ViHot.,**,. 

Lipschitz dependence on U and 0. 

Given 0 G i? 4 (r/ 4 ) and Ui, U 2 G x M), the difference between the solutions 

and Ve,d>,Ui fulhlls 

+ ^e,d>)'^{Ve,d>,Ui — Ve,<l>,U2) = (1 “ Xl){QE,<t>{X2U2 + 14 , 0 ,{ 72 ) ~ Qe,d>i.X2Ui + Ve^d>,Ui)) 

TNe,0(172) — + Pe,0(14,0,7/2) “ P£,0(14,0,t/i)• 

By dm), the terms involving satisfy 

||N£,0(17i) — Ne_,^(f/ 2 )||f; 0 ,a^jg 3 ^ < ce ^ ||172 — 

By dm), the terms involving can be estimated with the diherence between the 
solutions, that is 

I !?£, 0 ( 14 , 0 , 7 / 1 ) ~ ?£, 0 ( 14 , 0 , 7 / 2 ) I lc°’“(R3) — ce ^ ||I4,0,7/i ~ 14 , 0 , 7/2 1177^’“(R3), (83) 

and 


11 ( 1 - ~ Xl)(Q£, 0 (X 2 l 7 i + 14,0,7/1) ~ Qe,(l){X2U2 + 14,0,7/2)) I Ic“’“(R 3 ) — 
ce ^ ( 1 114,0,7/1 ~ 14,0,7/2! lc^'“(R3) + 11^1 ~ ^ 2 |lc'^.“(SexR))- 
Therefore, applying to the right-hand side of (l83|) . we obtain 

|| 14 ,0,7/1 ~ 14,0,7/2 II 4 
ce ^ (|| 14 ,0,7/1 ~ 14,0,7/2 I lc'^’“(M 3 ) + 11^1 ~ ^2! |( 74 ,a^j,^^jg^), 

thus, reabsorbing the norm of the diherence between the solutions, 

21114 , 0 , 7/1 ~ 14 , 0 , 7 / 2 1lc^’“(R3) 4 (1 ~ ce ^ )||14, 0 , 7/1 ~ 14 , 0 , 7 / 2 1Ic4.“(r3) 

< ce / ||f/i — I 72 I 

The Lipschitz dependence on 0 can be treated with a similar argument. It is worth to 
point out that also the potential Tg ^ depends on 4, through the approximate solution 
and the cutoh function. However, this dependence is mild enough for our purposes, in 
fact the diherence of the potentials Tg^^^ — Tg^^^ is exponentially small in e. 
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6.3 Invertibility in a neighbourhood of the linear problem 

Now we look for a solution to equation W!\ respecting the symmetries of the Torus. 
First we study the linear operator 

Proposition 20. LetQ<5< 0rX±l) and 0 e 54(1/4). For any f e x M), 

there exists a unique solution U = G^if) in x M) to C^U = f such that 

11^1 lc^’“(EeX]R) E C*! I/I 1(270. 

for some constant C > 0 which is independent of e. 

If / respects the symmetries of the Torus, then also the solution U = G^f does. In 
other words, G^ maps x M) into x M). This fact follows from uniqueness. 

It is useful to see that we can control the odd part of the solution with the odd part 
(in t) of / and the same is true for the even parts. 

Lemma 21. Let 0 < 6 < ^/W{1) and f e x M). Let U E x M) be the 

solution to C^U = f. Then 

I I l^ol lc/“(EeXR) E c| l/ol |(70 .q.(25x]R) 

[ I l^e| < c| |/e| |(270,a(.j2^^g^, 

where c is the constant found in Proposition W(\ 

Proof We set, for any {y,t) G x M, U{y,t) := U{y^-t) and f{y,t) := f{y,-t). 
Using that W is even and n* is odd, we can see that C^U = f. Therefore, subtracting 
and multiplying by 1/2, we get 

that is C^Uoiyjt) = fo. In addition, 

/ oo noo 

Uo{y, t)v[{t)dt = / fo{y, t)v^{t)dt = 0, 

-oo J —oo 

for any y E hence Uo = G^{fo), so in particular the hrst estimate holds true. The 
second one can be proved by a similar argument. □ 

Now we prove Proposition |20l with the aid of some Lemmas and Remarks. 

First we consider the spectral decomposition of We will denote by {^j,(f>j)j>o 
the eingendata of — A^. We observe that Aq = 0, Xj > Ai > 0, (fo is constant and, 
without loss of generality, we can assume that ||0j||L2(E) = 1 (see [26]). Similarly, 
we will denote by {p,k}k>o the eigenvalues of L* = —du + W"{v^{t)). In [2l], Miiller 
proved that fiQ = 0, and the corresponding eigenspace, that is the Kernel, is generated 
by v^(t), while /Xfc A > 0 (see also [20] )• 

Remark 22. The eigenvalues of are {p^k + £‘^Xj}j^k>o, thus al non-zero eigenvalues 
are positive and bounded away from 0, indeed pk + ^‘^Xj > e^Xi > 0. 
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Lemma 23. Let 


Ce : X R) ^ X 


be defined by the duality relation 


CsUu t /2 ) = / + dtU,dtU 2 + W (vfit))Ui [/2 \da{y)dt, 

/ JEeXR I J 


for any Ui,U 2 G C^’“(Se x R). Then 


Ker{Ce) = span{vfit)). 


Proof. It is possible to see that {Kj, 4>s,j)j>o ■= ^‘^4>j{^y))j>o ^-^e eigendata of 

and (fej are orthonormal in Any function w G x R) can be expanded 

in Fourier series as follows 


U{y,t) = ^Ufit)(t)e,j{y) 


where 


Uj{t)= U{y,t)(t),j{y)da{y). 


If C^w = 0, applying the operator to each term in the series, we get 
- duUfit) + Xs,jUfit) + w” {vfit))Ufit) = 0 
for any j > 0, so Uo{t) = cv'fit) and Wj = 0 for j > 1. 


0:= !^U e HfiE, xR): j 


U{y,t)vfit)da{y)dt = 0 


be the orthogonal to n*(t) in x 


Lemma 24. For any f G x R) satifying 


f{y, t)v'fit)dt = 0 for any y G 


there exists a unique U G x R) sueh that 


CeU = f 

IZc ^(y^ t)vfit)dt = 0 for any y G 
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Proof. At first we observe that 


l|t^ll= / \Vj:My,^)\" + iduU{y,t)f + W''{v:iz))U\y,z)da{y)dt (84) 

JEe xR 

is an equivalent norm on O, that is, for any 17 G X, we have 
C£,l||17||/fi(s^xR) < ||f^|| < C£,2||17||/fi(SjxR), 

for some constants Cs,i,Cs ,2 > 0. In fact, by the spectral decomposition of (see 
Remark |22]) , 


CJJUd(T{y)dt > e^Xi / U‘^da{y)dt. 


'SeXi 


' Se Xl 


Since hh”(n*(f)) is bounded, a pointwise estimate yields that 

f C^UUd(T{y)dt > f + {dttUYda{y)dt — c f U‘^da{y)dt, (85) 

7SeXR 7SeXR JScXR 


for some constant c > 0. Now we point out that, for any 0 < A < 1, we have 


C,UUda{y)dt = X / C,UUda{y)dt + {1 - X) 


C^UUda{y)dt > 


'SeXR 


'SeXR 


x( [ + {dttUyda{y)dt-c f U^da{y)dt] + {1 - X)e^Xi f U^da{y)dt, 


V^SeXR ^SjXR / ^EgXR 

SO, in order to prove the lower bound, it is enough to choose A < £^Ai/(c + £^Ai). As a 
consequence, by the Riesz representation theorem, for any / G -^^(Eg x M) such that 


f{y,t)v*{t)da{y)dt = 0, 


( 86 ) 


the equation Cell = f admits a unique solution U G O. We observe that orthogonality 
condition (l86|) is necessary for solvability, since 



f{y,t)v'^{t)da{y)dt 


/ CeU{y,t)v^{t)da{y)dt = 

'SeXR 

U{y,t)Cev[{t)da{y)dt = 0. 


If in particular / satishes flM)) . then, by proposition 8,4 of [2S], also w satishes fIM)) . □ 

Now we are ready to conclude the proof of Proposition l20l 

Proof. There are two more steps. As hrst we need some regularity theory to estimate 
the (7^’“ (Eg X R) norm of the solution t/ if / G 8^’°‘{'Ee x R), then we have to iterate 
the estimates to deal with the operator C^. For the hrst step, see Proposition 8,3 of 
|26) . As regards the second one, we argue as follows. 
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If / G X M), the above discussion yields that we can hnd U G x M) 

such that 


\^eU = f 

\ll^llc^’“(SeXR) — ^ll/llc°’“(SeXR)> 

for some constant C > 0 independent of e. Now, by the same argument, we can hnd 
U G X M) satisfying 

UeU = U 

\ 11^1 lc^’“(SeXR) — ^11^1 lc°’“(SeXR) — ^1 I/I lc“’“(EeXR)’ 

for some constant C > 0 independent of £. To conclude the proof, we have to show 
that U G X M) and 


Cj’“(EeXR) — ^1 I/I lc°’“(Ee. 


(87) 


In order to do so we apply a bootstrap argument. We differentiate fl87|) with respect 
to Hj and we get 


= u,. 


By fl87|) . we get that Uj G x M) and 

I IT/'I lc2’“(EeXR) — ^1 l^il lc°’“(EeXR) — ^11^1 lc^’“(EeXR) — ^1 I/I lc'°’“(Ee> 
In the same way, taking the derivative with respect to t, we get 


C,Ut = Ut--W {v,{t))v,{t)U. 


Exactly as before, we have 


11^*1 lc^’“(EsxR) — ^(11^/ lc°’“(EsxR) 11^ ('*^*(^))W(^)77| |c'0.“(S 

^(11^1 lc5’“(EeXR) + 11^1 lc°’“(EsXR)) — ^(1 I/I Ic5’“(EsXR) + 11^1 lc^’“(EsXR)) — ^ll/llc°’“(E 


xR)) — 


Therefore we have 


\\v%u^ps)\\o. < < c||/||^o..(^, 

Differentiating the equation once again, we get 

||V^(7/^,)||oo + [V\UiJs)]a < C^II/llcr(E.xR)- 

In conclusion, we have fl87p . 


□ 
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6.4 The proof of Proposition [TOt solving equation ( l50l) by a 
fixed point argument 

Equation flSU]) is equivalent to the fixed point problem 

U = T,{U) := G,| - +p(|/)n:(t)|. 

Once again, we will solve it by showing that T 2 is a contraction on the ball 
A 2 := {U € X R) : I|(7|Ic,*..,!:.x«) < 

provided 02 > 0 is large enough. First we observe that, by definition of p, the right 
hand side is orthogonal to for any p G thus we can actually apply the operator 
Gg. Moreover, if U respects the symmetries of the Torus, then also the right-hand side 
does, thus, applying G^, we get once again something that respects these symmetries. 
Now we show that, if I|G|< 02 ^^, then also T 2 {U) satisfies the same upper 
bound, for some large constant G 2 . 

We note that 


||X4fo('?^£,</))|lcO'“(SgxK) — ^ ’ 

for some constant c depending just on fT, r and the geometric quantities of S, and the 
same is true for p{y)v^{t). The other terms are smaller, for instance, using (l3^ and 
the fact that V is exponentially small, 

||XlQ£, 0 (fo + fo)llcO'“(SeXR) ■ 

Similarly, we can see that In addition, since all the 

coefficients of are at least of order e, we get that 

||X4R£,(/>(G)|Ic<0.“(S5xR) — ^11^1 lc^’“(SeXR) — ^ ■ 

As regards the Lipschitz dependence on U, we observe that 


\\xi{QeAUi + v) 


QeAU2 + vm 




<C£='||Gi-G2 | 


CAi^e 


and 


I |X4(R'£,</)(Gi) ^£, 0 (^ 2 )) I IcO.“(SjxR) — *^^11^1 ^ 2 ! Ic^.a^^jxR)- 

Estimate of the odd part of the solution 

Up to now we have proved the existence of a solution Ue,(j, to equation fl5n|) sat¬ 
isfying I |G£, 0 | < ce^. However, we point out that the only terms of or¬ 

der in the right-hand side come from X4F{veA- fact, as we observed above, 
||T(G, 0)|lc°’“(SexR) — particular the same is true for 
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Moreover, by Proposition [151 

/ oo poo poo 

XiF{ve,^)iy,t)%{t)dt= / F{v,^^){y,t)v[{t)dt+ / {x4 - l)F{v,^^){y,t)v[{t)dt 

-OO J —oo J —oo 

is of order since the second term is exponentially small. Going back to Section 5, 
it is possible to see that the only terms of order in are even in f, thus the 

odd part of the right-hand side is of order and therefore, by Lemma 1211 the same 
is true for namely \\Ue,<p\\cf^^Y:exR) ^ 

Lipschitz dependence on (p. 

Let us £x 01,02 € A 2 . To simplify the notation, we set, for fc = 1,2, 14 := 

Uk := Ue^ 4 ,^, 14 := so on. In this proof, £ will always be small but fixed, 

and we will be interested in the dependence on 0. 

First we note that 


X 4 (T(hi) - F(h 2 )) = T(hi) - F{v,) + (X4 - l)(T(hi) - F{v 2 )) 
The first term satisfies 


||F(hi) -F(h2)|lc0“(SexR) < c£3|0i -02|c4.-(s), 
because, for instance, 

k^(|Vs0l|^ — |Vs02n'yi^^| < C£^(|Vs0l| + |Vs02|)|01 — 02|c4.“(S) < C£^|01 “ 02|c4.“(S) 

The other terms are similar, or even easier to treat because there is already an that 
multiplies everything (see section 5,1). The Lipschitz constant of the second term is 
exponentially small in e, namely 


||(X 4 - l)(T(hi) - T(h 2 ))|lc 0 “(s,xR) < ce “^401 - 02 |c 4 ,-(s)- 

Using the Lipschitz dependence of V on the data proved in Proposition [H] and the 
dehnitions of 0 , and Re, 0 , it is possible to see that 


||Mi(Ui) - M2(U2)|lc0“(SeXR) - ~ ^2|lc-0“(SeXR) + ~ ‘/’S | C4.“(S)), 


||X 4 (Qi(Ui -|- Vi) — <52(^2 + 14 


IC^’“(SeXl 


<ce^{\\U4-U: 


2llc^’“(SeXl 


+ 101 ~ 02|c4.“(S)), 


||X4(Ri(Ui) - R-2(U2))||c4.<^(s^xR) < C£||Ui - U 2 | + €£"^101 - 02|c4.“(S)- 


Now it remains to deal with p{y), that also depends on £ and 0. We write, for any 
y eT,^ and 0 e S4(r/4), 


/ OO 

F{ve,^){y,t)v[{t)dt + {pi( 0 )(l/) +P 2 ( 0 )( 2 /) +P3(0)(2/) +P4(0)(2/)}wl(f), 

•OO 
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where we have set 


1 /‘°° ! 

Pi( 0 )(|/) := — / (1 - Xi)F{ve^^){y,t)v^{t)dt, ( 88 ) 

J —OO 

1 Z*^ / 

P 2 i<P)iy)-=— xiQeAU+ V){y,t)vl{t)dt, (89) 

.7 —OO 

1 / 

P3i<P)iy)-=— XiM£,^(t^)(|/,t)i;l(t)rft (90) 

.7 —OO 

1 Z*^ / 

P4(0)(2/) := — / X4R£,<A(t^£,0)(2/,^)wl(^)c?^ (91) 

J—OO 


and U := V ■= Since we want to deal with fnnctions dehned on S, we will 

set, for any y G pj(0)(£|/) := Pi(0)(|/), for i = 1,..., 4. It follows from Proposition 
m and that 


F{ve,^J{y,t)v'^{t)dt - / F{ve,4,^){y,t)v[{t)dt 


< C€^\(t)i - 02 I 


(S)- 


(92) 


In addition, by the previons discnssion. 


bi(0)|co.“(s) < ce 

|Pi(0i) -Pi(02)|co.“(s) < ce““/00i -02 |c4,-(S)- 


(93) 


Fnrthermore, by the Lipschitz dependence of V on the data, proved in Proposition |9l 
and by the fact that I |f^| < C 2 £^, we have 


I b2(0)|co.“(s) < 

[|P2(0l) - P2(02)|cO.“(S) < C£^(|01 - 02 |c 4.-(S) + llh^l “ ^^2 11 c0“ (S^ xR))' 

and, similarly 


(94) 


||P3(0)|co.-(s) < ce 

\|P3(0l) - P3(02)|cO’“(S) < Ce““/^(|0i - 02 |c 4.“(S) + H^l “ ^^211 c0“(Se xR))' 

As regards p 4 , we give a hrst, rongh estimate that is enongh to prove the Lipschitz 
dependence of 1/ on 0. However, we will see later that this estimate is actnally not 
enongh to solve the bifnrcation eqnation, thns we will improve it in Lemma [25], nsing 
the estimate of the odd part of f/ (see section 7). 


|p4(0)|co.a(s) < 

|P4(0l) -P4(02)|cO.“(S) < ce:||[/i - fl2||c0“(SeXR)- 


(96) 


In conclnsion, the eqnation satished by the difference of the solntions Ui — U 2 is of the 
form 


4(^1 - U 2 ) = gi(j)i)iy,t) - g{(j)2){y,t), 
where g{(j)i) and Ui satisfy 

POO POO 

/ ig{<Pi) - g{h))iy,t)%it)dt= {Ui-U 2 )iy,t)vl{t)dt = o, 
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thus, by Proposition 1201 

||h^l - ^2||c4.“(EeXR) ^ Ce\\Ui - t/2||c4.-(s,xR) + “ 02|c4.-(S), 

and hence, reabsorbing the hrst term of the right-hand side, 

- t^2||c'4.“(EeXR) ^ - 02|C4.-(S)- 


7 Solving the bifurcation equation 

7.1 The proof of Proposition [TT1 

First let us hx some notation. For any 0 G and 0 < e < 1, 0 I 3 will be 

the volume of the interior of that is its 3-Lebesgue measure. Moreover, we set 

Bi := {x = Zs{y, t + (t){ey)) : -5 - r < f < 0} 

B 2 := {x = Z^{y, 2 ;) : 0 < f < 5 -|- r/2e}, 

Vi will be the volume of B^, for i = 1, 2, and A := M.^\B. Now we note that 


{1 - Ve,^{x))dx = / {1 - v^^^{x))dx + / {1 - v^^^{x))dx 
* Ja Jb 


and 


/ (1 — v ^^^{ x))dx + / Idx — 2 (| Ee_^|3 — Vi ) -f Vi -|- V2 — 2 | Ee^|3 -f V2 — Vi . 

JA JB 

In the forthcoming integrals, we will use the natural change of variables induced on 
Vt/i; by the parametrization ^^(y) = e~^Y{ey) (see flM|) h The absolute value of the 
Jacobian determinant is e‘^\^{z + cos(£yi) -1- {z + e~^)e~^\/2'^, thus we can see that 


1‘2-k je 


|S£,,^|3 = 27re M dy^ {{t + (p{ey^) + e ^)^cosyi 

Jo J-l/e-<t>{eyi) 

+{t + (p{ey^) + e~^)e~^V2}dz = e~^2TT^V2 + e~‘^ f (p{C)da{C) 


(97) 


r-27r 


-|-27re ^ / 0^('d)(cos'd-I-\/2/2)(i'd-I- 


27r 


f27r 


0^('d) cos{'d)d'd, 


since the surface integral over of a function -0 of the variable y^ is given by 

(‘2w/6 


I ^{y)d(T{y) = 27re M (cos(£yi) + V2)^/J{y^)dyl. 
'Se Jo 


(98) 
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Similarly, we can show that 


27re: 


-1 


^‘I'Kje pQ-\-Tl2e 

dji / 


r-27r/e ^0 


1^2 - = 
£2{(t + 0(£yi) +£“^)^cos(£yi) 
+{t + (j){eyi) + e~^)e~^^/2]dt 


(99) 


- 27 r£ W dy^ e'^{{t + (j){ey^) + e ^)^cos(£yi) 

Jo J-6-r/2e 

+ (t + 0(£yi) + £“^)e:“^\/2}c?t = 

pQ-^Tl2£ n2'K 

27re“^ / tdt / {2\/2 + 4e0('i9i) cos('(9i)}(i'(9i. 

Observing that 

= 'Ce,<^(£yi, t) + (1 - X5(a^))(H(x) - h£,0(£yi, t)) 

we compnte 


' B 


( 100 ) 


n£,0(x)ci^lOl) 


f27r/e 


f 6 +r/ 2 £ 


.-1 


= 27re W rfyi / t){(t + ^(eyi) + £ ^)^cos(eyi) 

Jo J-O-rj^e 

+{t + 0 (£yi) + £“^)£“^\/ 2 }rft 

r*27r/£ p6+rl2£ 

+2Tie~^ I dy^ (1 - X5)(H(x) - ^^^^(eyi, t))£^{(t + ^(eyi) + £“^)^ cos(£yi) 

Jo J-6-Tl2e 

+(t + 0(£yi) + e~^)e~^V2]dt. 

The second integral is exponentially decreasing in e, and the same is trne for its Lips- 
chitz constant. As regards the second one, exploiting the symmetry of n*, r] and of the 
domain, we can see that 

i'2TT/e i'6+tI2s 

2Tre~^ dyi / e^v,^^{ey^,t){{t + (j){ey^) + e~^Y cos{ey^) 

Jo J-6-Tl2e 


+ {t + (j){ey^) + e ^)e: ^\/2]dt 


2718 


-1 


f6+r/2£ 


f27r 


tVi,{t)dt / {4£0('di) cos('d) + 2 \/ 2 }(i'di + (^£( 0 ). 


with Gl satisfying fl55|) . Thns, taking the snm of fl98|) . ([99]), fllOip and fll02p . 

/ (1 — n£,</,(x))(ix = £“^47r^\/2 + 2 e“^ / 4>{C)d( ( 102 ) 


+27re 


-1 


i‘6+t/2s 


-27r 


t{l - Vi,{t))dt / {4£0('di) cos('di) + 2\/2}d-di + Gg(0), 


with Gl satisfying fl55|) . 
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It remains to deal with the term involving 

/ \We,<f>{x)\dx = / \we, 4 >ix)\ipsix)ip-s{x)dx 

JR3 Jr3 

< c||we,0|L4.a.^3s / (p_s{x)dx < ce^ 

^ ^ Vr3 

and, by Propositions |9] and [TOl 

/ {We,cl,:^ix) - We,<l,2i^))dx 
Jr3 

^ ^£,021 Ic'^’“('r3') / ^—s^x'jdx ^ c£ 101 02|c^’“(s)) 

^ ^ Vr3 

for any 0 i ,02 e C‘^’“(S )3 satisfying |0i|c4,-(s), |02|c4.“(s) < c£. 


7.2 The proof of Proposition \n\ 


Before giving the proof, we state a technical Lemma, in which we prove that the 
term is small enough. 

Lemma 25. For any £ > 0 small enough, for any c > 0 and for any 0, 0i, 02 satisfying 
|0|c4.“(s)5 |0i|c4.“(s), |02|c4.“(s) < ce, we have 

\ |P4(0)|co.“(s) < ce^ 

\|P4(0l) — P4(02)|cO'“(S) < C£^|0i — 02 |cO’“(S)5 
for some constant c > 0. 

Proof We write = (t4,0)o+(t4,,^)e. By Proposition HOl we know that 11(^4,</.)o|lc0“(SexR) ^ 
c£^, therefore 11-Re,0((f4, 0 ) 0 )|lc'4-“(SexR) 4 ce^, since all the coefficients of are at 
least of order £. It remains to deal with the even part Ue- We will see that all the 
terms of order in the expression of R£,0(f4,0)e will vanish after projection. This can 
be seen by a direct computation 

X4ReAiUe,^)e) = SxAH^'" (f/,,0)et) 

+Hdut{U,Ae + W''{vAHdt{U,Ae + a\^dAUeAet) 

+ iHdtiU,Ae + aYdAUeAet)RsiUeAe + 4,0((R..0)e) } , 

where ReAiR^Ae) some linear operator with coefficients of order at least All 
the terms of order e are odd, thus they vanish when we multiply by and integrate, 
the other ones give rise to terms of order e^, being 17^,^ of order e^. □ 

Now we are ready to prove Proposition [T2l 


Proof. In view of Proposition [151 the system of equations fl52|) and fl55|) is equivalent 
to the hxed point problem 


0 = T3(0) :=-p(^C + £ ^{pi(0) + P2(0) + P3(0) + P4(0)}, 


4\/27r^e 


f6+T/2e 


t(l — vA))dt + e^Ge{ 4 >) 
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where P : x M ^ is the projection onto the hrst component. We 

will show that T 3 is a contraction on the ball 

A3:={0eC'"’“(S), :|0 |c4,.(s)<C'3£}, 
provided C 3 is large enough. 

Using once again the same estimates as in the proof of Proposition [10] and the fact 
that Lipschitzianity of U with respect to 0, we can see that pi and are exponentially 
small in e, that is they satisfy, for instance 

f |Pi(0)lco.“(s) < ce-“/^ 

\|Pl(0l) -Pl(02)|cO.“(S) < Ce-^l^\px -02|c4.“(S), 
for any 0, 0i, 02 £ A 3 . Similarly, by (l39|) . we can see that 


b2(0)|cO.“(S) < C£® 

|P2(0l) — P2(02)|cO’“(S) < C£®|01 — 02 |c4.“(E)- 


The term is small according to Proposition IT51 The most difficult term is the 

one involving since there are some coefficients of order e: and P is just of order e^. 

However, we verihed in Lemma [25] that these terms do not give rise to terms of order 
after projection, thanks to the symmetries. 

The second component can be treated in a similar way. In fact 


4\/2' 




f6+T/2£ 


f(l — Vi,{t))dt < ce 


and it is independent of 0. To conclude, e^Gs^cf, is small according to Proposition [TT] 
In conclusion, T 3 is a contraction of the ball A 3 , provided C 3 is large enough. □ 
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